A numerical study of one-patch colloidal particles: from square-well to Janus. 
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We perform numerical simulations of a simple model of one-patch colloidal particles to investigate: 
(i) the behavior of the gas-liquid phase diagram on moving from a spherical attractive potential to 
a Janus potential and (ii) the collective structure of a system of Janus particles. We show that, for 
the case where one of the two hemispheres is attractive and one is repulsive, the system organizes 
into a dispersion of orient at ional ordered micelles and vesicles and, at low T, the system can be 
approximated as a fluid of such clusters, interacting essentially via excluded volume. The stability 
of this cluster phase generates a very peculiar shape of the gas and liquid coexisting densities, with 
a gas coexistence density which increases on cooling, approaching the liquid coexistence density at 
very low T. 
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I. INTRODUCTION 

The synthesis of colloidal particles with controlled 
anisotropy is central in today research. One of the main 
ideas is the development of a set of colloidal molecules [1 - 
[7] which can be useful to generate, on the nano and 
micron-scale, collective phenomena presently observed 
only at atomic or molecular scale as well as additional 
phenomena, induced by the possibility of controlling the 
parameters of the effective interaction potential between 
colloids. The self-assembly of a colloidal diamond crys- 
tal, to be used in photonic applications [8 , is one of these 
technological relevant goals. The anisotropy can be in- 
duced not only by building colloidal molecules (i.e. col- 
loids with peculiar non-spherical shapes) , but also (with a 
promising alternative) via the process of patterning the 
particle surface, generating in this way particles inter- 
acting in very different way according to their relative 
orientation[Tll9Hl3]. 

This vision of a material science, based on the design 
and self-assembly of materials with required properties 
has stimulated, beside the experimental work, a large 
amount of theoretical and numerical studies based on 
simple, primitive, anisotropic potentials [14- 20 . Indeed, 
the possibilities of modifying the surface properties are 
endless, including the number of patches, their width, 
their location, their chemical specificity, and seem to pose 
no limits to the design of specific particle-particle inter- 
action potentials. In this large parameter space, primi- 
tive potentials — modeling repulsion as a hard-core and 
attraction as a square-well interaction — can provide 
a useful reference system to deeply investigate the role 
played by the number of patches, their width, their spa- 
tial location and the role of the attraction range. In addi- 
tion, there is some consensus that these studies can shed 
light on the aggregation properties of proteins and on 
the sensitivity of the aggregates on the protein surface 
properties [2TH24] . 

Another useful aspect of the study of primitive mod- 
els is the possibility to accurately investigate theoreti- 



cally and numerically their phase diagram, casting the 
self-assembly process of these systems into a wider ther- 
modynamic perspective. Important questions about the 
nature of the self-assembly process and its reversibility, 
the relative stability of aggregates of different sizes, the 
competition between self-assembly and phase-separation 
can in principle be addressed with an accurate study of 
these models. Finally, in some cases, numerical simu- 
lations can be compared with state of the art integral 
equation approaches for non-spherical potentials provid- 
ing a benchmark for the possibility of developing a fast 
and accurate prediction of the structural properties. Sev- 
eral efforts in this directions, capitalizing on studies of 
molecular associations [25- 28 , are taking place in these 
days[: 



In this article we report a study of the phase behav- 
ior of a very simple primitive potential, proposed in 2003 
by Kern and Frenkel[14 with the aim of exploring the 
scaling of the critical parameters (including the reduced 
second virial coefficient) on the number and width of the 
patches, and the possible implications for the phase be- 
havior of globular proteins. We limit ourselves to the case 
where the surface of the colloid particle is divided into 
two parts, respectively repulsive and attractive. On de- 
creasing the surface area corresponding to the attractive 
part, the potential interpolates between the well-known 
isotropic square well potential and the symmetric case of 
an evenly divided surface, commonly indicated as Janus 
potential. We find that an unconventional phase dia- 
gram characterizes Janus particles, due to the onset of 
a micelle formation process which takes place in the gas 
phase, providing additional stability to this phase as com- 
pared to the liquid one. We characterize the structural 
and connectivity properties of the system in a wide range 
of temperatures T and number densities p, to clarify the 
origin of the aggregation process and the mechanisms be- 
hind the stability of the cluster phase. Our study demon- 
strates how a small change in the attractive surface has 
profound consequences on the collective behavior of the 
system. 



II. MODEL AND SIMULATION TECHNIQUES 

The Kern-Frenkel potential is a paradigmatic model 
for highly anisotropic interactions. In this model, a hard- 
sphere of diameter a is complemented by a set of unit vec- 
tors {rii}, locating the position of the center of a patch on 
the particle surface. Each patch can be reckoned as the 
intersection of the sphere with a cone of semi-amplitude 
and vertex at the center of the sphere. In the case stud- 
ied in this model, each particle has only one patch. When 
unit vectors the patch fii of particle i and fij of particle 
j form an angle smaller than #, and in addition the dis- 
tance between the center of the two particles is between 
cr and a + A, then an attractive interaction of intensity 
uq is present. More precisely, the two body potential is 
defined as: 



u(r ij )=u sw (r ij )f(Q ij ) 



(1) 



where u sw (rij) is an isotropic square well term of depth 
uq and attractive range a + A and f(&ij) is a function 
that depends on the orientation of the two interacting 
particles Qij. The angular function f(£lij) is defined as 
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(2) 
The diameter of the particles and the depth of the 
square well have been chosen as units of length and en- 
ergy respectively, i.e. a = 1 and uq = 1. 

In practice, two particles interact attractively if, when 
they are within the attractive distance <r+ A, two patches 
are properly facing each other. When this is the case, 
the two particles are considered bonded. The fraction of 
surface covered by the attractive patches x 1S related to 
by the relation x = 1 ~ c 2 os ^ 

Structural properties of the system have been evalu- 
ated by mean of standard MC simulations, for a system 
of N = 5000 particles. Extremely long simulations, of the 
order of 10 9 MC sweeps, have been performed to reach a 
proper equilibrium state of the system. Here a MC sweep 
is defined as an attempted random translation and rota- 
tion for each particle. To calculate the location of the gas- 
liquid critical point we perform grand canonical Monte 
Carlo (GCMC) simulations [32], complemented with his- 
togram reweighting techniques to match the distribution 
of the order parameter p — se with the known functional 
dependence expected at the Ising universality class crit- 
ical point [33]. Here e is the potential energy density, p 
the number density and s is the mixing field parameter. 
We did not performed a finite size study, since we are 
only interested in the trends with \. We have studied 
systems of different sizes, up to L = 15. For each studied 



X — using the methods described in [34] — we calculated 
the critical temperature T c and density p c for values of 
cos# between —1 and (at fixed A = 0.5). Temperature 
is measured in reduced units, i.e. Boltzmann constant 
k B = I- 

We also performed Gibbs ensemble simulations to eval- 
uate the coexistence curve. The GEMC method was de- 
signed [35 to study coexistence in the region where the 
gas-liquid free-energy barrier is sufficiently high to avoid 
crossing between the two phases. Since nowadays this 
is a standard method in computational physics, we do 
not discuss it here. We have studied a system of (total) 
350 particles which partition themselves into two boxes 
whose total volume is 2868a 3 , corresponding to an av- 
erage density of p = 0.122. At the lowest T this cor- 
responds to roughly 320 particles in the liquid box (of 
side « 8a) and about 30 particles in the gas box (of side 
« 13a). Equilibration at the lowest reported T required 
about one month of computer time. 

The model, in the case x — 0-5, can be related 
to the experimental system that is currently under 
investigation [TT] . In these newly synthesized Janus parti- 
cles, the repulsive interaction has an electrostatic origin 
and the attractive part is hydrophobic At the present 
time, experiments by Granick's group [11 have focused 
on the analysis of the structure of the aggregates sedi- 
mented on the bottom surface due to gravity. Micelle 
formation has been observed. While the interaction range 
in the Granick's experimental system is about 0.1 of the 
particle size, nano-sized particles synthesized with the 
same protocol would indeed give rise to potential ranges 
similar to the one we have selected. It is also in prin- 
ciple possible to modify the range of the interaction by 
tuning the physical properties of the solvent (for exam- 
ple the ionic strength or the solvent dielectric constant). 
Another interesting possible experimental realization of 
longer range interactions can be achieved by mean of the 
recently measured Casimir critical forces [for a highlight- 
ing review see for example Ref. [36 or Ref. [37]] dissolving 
the particles close to the critical point of the solvent. By 
tuning the distance from the critical point, the range of 
the interaction can be controlled and tuned close to the 
value we have explored. Moreover, in the critical Casimir 
effect, the coating of the two hemispheres will controls 
if the interaction is attractive or repulsive. Other pos- 
sibilities of probing different ranges will be offered by 
studies of Janus-like proteins. Indeed, the hydrophobin 
proteins (extracted from fungi) are good candidates for 
nano- Janus particles, as discussed in Ref. [38 J. Even in 
this protein case, the bulk phase diagram has never been 
explored. 



III. RESULTS: FROM SQUARE- WELL TO 

JANUS 



The Kern-Frenkel model offers the possibility to con- 
tinuously change the coverage interpolating from the 
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FIG. 1: Phase diagram of the one-patch Kern-Frenkel po- 
tential with attractive range 0.5a for different values of the 
coverage x> interpolating between the square- well potential 
(X — 1) an d the Janus potential (x = 0.5). 



isotropic square-well potential to the symmetric Janus- 
like one, when the coverage moves from ^ = ltox = 0.5. 
To investigate how the phase diagram of Janus particles 
arises, we start by looking how the gas-liquid coexistence 
is modified on progressively reducing x- Fig. [T] depicts 
the gas-liquid phase coexistence for several x values, ex- 
tending the original data by Kern and Frenkel[T4]. While 
the critical density presents a significant decrease with 
decreasing x, the density of the liquid branch does not 
show any significant reduction, consistent with the pos- 
sibility of forming six or more bonds with neighboring 
particles when x > 0.5. Recent studies on the role of 
the valence (defined as the maximum possible number of 
bonded nearest neighbors) [I5j [39] , have suggested a pro- 
gressive reduction of the critical temperature T c and crit- 
ical density p c on decreasing the valence. Indeed, when 
the valence decreases below six, a significant reduction 
of the density of the liquid branch has been reported [16]. 
Estimates of the critical parameters have been calculated 
using grand-canonical simulations, i.e. simulations in 
which the chemical potential /i, T and the volume V are 
kept constant. In GCMC simulations, the number of par- 
ticles fluctuates. When p and T are close to their critical 
values, the number of particles fluctuates widely (since 
the compressibility, which is a measure of the variance 
of the density fluctuations diverges). The distribution 
of sampled densities, close to the critical point, follows 
an universal curve which depends only on the class of 
universality of the critical phenomenon. 

The critical parameters, resulting from the grand- 
canonical simulations, are summarized in Table [T[ The 
last column indicates the size of the largest simulation 
box, which has been progressively increased to compen- 
sate the decrease of the critical density and the associated 
shift towards smaller number of particles of the density 
fluctuations. Table [I] also reports the critical chemical 



FIG. 2: Coverage x dependence of the critical density p c and 
temperature T c . 



potential and the value of the second virial potential B%, 
normalized to the hard-sphere value B^ 
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also evaluated at the critical point. Previous work [19] has 
shown that B^/B^ becomes smaller and smaller with 
decreasing valence. The actual value of B^/B^ can 
provide an estimate of the effective valence of the system. 
A comparison with Table 1 of Ref . [19] suggests that when 
X reaches the value 0.6, the effective valence becomes less 
than 4 and that when x = 0-5 the effective valence is 
still larger than three. The critical parameters T c and p c 
are also shown graphically in Fig. [2] In agreement with 
the progressive reduction of the valence, T c and p c do 
decrease with decreasing \. Note that results presented 
in Figure [I] have a counterpart in the case where the 
attractive part is spread over two patches distributed at 
the opposite poles of the sphere (see Fig. 2 in Ref. [3T] ) . 



TABLE I: Critical parameters as a function of the coverage. 
The last column indicates the size of the largest studied sim- 
ulation box. 
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When x — 0.5 (the Janus case), a new interesting 
phenomenon appears. Fig. [I] shows indeed that, differ- 
ently from the standard behavior (and for all the studied 
X > 0.5 cases) the density of the gas phase along the 
coexistence line increases progressively on cooling. This 



peculiar behavior is discussed in full details in the follow- 
ing sections. For smaller values of x we have not presently 
been able to evaluate the phase diagram for two differ- 
ent reasons: i) the temperature region where the critical 
point is expected (around T = 0.17 for x — 0-4, by a 
quadratic extrapolation of the data reported in Table [I]) 
requires significant computational resources, ii) prelimi- 
nary tests have detected the formation of lamellar phases 
for x — 0-4 already at T > 0.17. Despite the disappear- 
ance of the liquid phase as an equilibrium phase in the 
presence of anisotropic potentials is a potentially rele- 
vant issue [40 , we can not at the present time address 
this point for the present model. A study of the stabil- 
ity of the liquid phase as compared to the (unknown) 
ordered phases will be the argument of a future study, 
which will require the use of algorithms to identify the 
possible crystal structures [4TJ [42] as well as free-energy 
evaluations to establish the relative stability compared to 
the liquid phase. A similar study for tetrahedral patchy 
particles has been recently reported [40, 43 . Indeed, it 
could well be, as suggested in a recent study of particles 
with two patches [31], that the reduction of the bond- 
ing angle could play a role analog to the reduction of 
the range in spherically interacting potential [44, 45 and 
limited valence potentials [40], where the liquid phase dis- 
appears when the range is smaller than w 0.2<r. 



IV. RESULTS: JANUS 

The gas-liquid coexistence for the Janus (x = 0.5) par- 
ticles is enlarged in Fig. [4] As discussed in a preliminary 
publication [46 , the phase diagram has a very odd behav- 
ior. The gas coexisting density, which typically is a de- 
creasing function of the temperature, here increases pro- 
gressively on cooling, approaching the coexisting liquid 
density. In a simple liquid, coexistence between gas and 
liquid is established on the basis of a compensation be- 
tween the gas and liquid free energies. The lower energy 
of the liquid phase is compensated by a larger entropy of 
the gas phase, which is acquired by significantly increas- 
ing the volume per particle. As discussed in Ref. [46] (and 
detailed more in the following), the uncommon behavior 
observed in the Janus case arises by a completely dif- 
ferent compensation mechanism between the liquid and 
the gas. The gas becomes the energetically stable phase 
(due to the formation of orientationally ordered aggre- 
gates, micelles and vesicles) and the liquid phase instead 
is stabilized by the larger orientational entropy of the 
particles. This odd behavior give rises to a gas-liquid co- 
existence curve in the P — T plane which is negatively 
sloped and to an expansion of the system on crossing 
from the gas to the liquid phase on cooling along isobars 
(see Fig. 4 in Ref. [46]). 

This anomalous thermodynamic behavior arises from 
the progressive establishment in the gas phase of clusters 
of particles which — due to the surface pattern proper- 
ties of the particle — organize themselves in particularly 



stable structures. Typical cluster shapes for different val- 
ues of the cluster size s are shown in Fig. [3] For small 
cluster sizes (s < 20) clusters are of micellar type, i.e. 
formed by aggregates in which the attractive part con- 
stitutes the core of the aggregate. For larger size, the 
cluster organization changes in favor of a double layer 
structure, reminiscent of vesicles, in which the inner and 
outer surfaces are repulsive and the inner core is attrac- 
tive. Here, and in the rest of the manuscript, clusters 
are defined as set of particles connected by an uninter- 
rupted path of bonds, where we define bonded any pair 
of particles whose pair potential energy is — u$. 

In this article we explore in details the properties of the 
clusters which develop in the gas phase, their structure, 
energy and abundance. We then investigate the gas and 
liquid phases with the aim of characterizing the collective 
structure of the system (both in real and wave-vector 
space) as well as its connectivity properties (percolation) . 



Critical Micelle Concentration 



The onset of micelles can be demonstrated and quan- 
tified by the study of the relation between the density 
pi of particles in monomeric state (i.e. un-bonded par- 
ticles) and the system density p. Fig. [5] shows that for 
T < 0.28, a sharp kink separates the ideal gas behavior 
(where all the particles are in a very dilute monomeric 
state and p\ = p) from a rather insensitive density de- 
pendence of the number of monomers in solutions, a clear 
indication that, at low T, the addition of particles to a 
constant volume system promotes the formation of ad- 
ditional aggregates. This behavior is indeed typical of 
micelle forming systems [47], and the location of the kink 
provides an estimate of the critical micelle concentration 
(am. a), which in the present case varies from p = 10 -3 
down to p = 10 -4 when T changes from T = 0.28 to 
T = 0.25. 

One may wander why this system, on cooling, does not 
show a typical gas-liquid coexistence. In simple fluids, at 
low T, the gas phase does not show significant clustering. 
Indeed, on increasing the density, the system transforms 
into a liquid phase, by establishing an infinite size perco- 
lating clusters, in the attempt to minimize the potential 
energy, being the entropic loss in the free energy (asso- 
ciated to the restricted sampling of the system available 
volume) made less relevant by the small T. In the present 
case, the Janus potential allows for the establishment of 
significantly bonded aggregates, i.e. with a small poten- 
tial energy, without the need of forming an infinite size 
cluster. By exposing the hard-core part to the exteriors, 
these clusters do not feel any driving force toward further 
clustering. 
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FIG. 3: Typical cluster shapes of different size, extracted from simulations at T = 0.27. 
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FIG. 4: Phase diagram of the Janus particles. Filled (red) cir- 
cles indicate the gas-liquid coexistence lines with the (blue) 
triangle denoting the critical point. The filled and open 
squares indicate the percolating and non-percolating state 
points, respectively, whereas (green) diamonds indicate the 
simulations that show a lamellar phase. Dashed lines connect 
coexisting state points whereas the two dotted lines refer to 
the two paths followed in the calculation of the structure fac- 
tor of Fig. [8] 



B. Cluster Size Distributions, percolation 

To quantify the effect of clustering and its p and T 
dependence we show in Fig.|6|the cluster size distribution, 
N(s). 

Below the critical temperature, the gas phase becomes 
populated by rather stable clusters. Particularly signif- 
icant are clusters of size between 10-15 (micelles) and 
clusters between size 40 and 50 (vesicles). Other sizes 
are also found, but their statistical relevance decreases on 
decreasing T. This feature appears clearly in the cluster 
size distribution, N(s), shown in Fig. [6] On decreas- 
ing T, the monotonic decaying N(s) develops a shoulder 
around T c , which evolves into a clear two-peaked function 
on further cooling, signaling the appearance of micelles 
and vesicles (Fig. [6| (a)). Part (b) of Fig. [6] shows the 
evolution of N(s) on increasing p. The micelle peak pro- 
gressively empties in favor of the vesicles peak, which at 
low T and large p are the most stable structures. Similar 
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FIG. 5: Relation between the monomer number density pi 
and the system number density p at different T. The flat re- 
gion for T < 0.28 indicates that the addition of monomers 
preferentially results in the formation of aggregates. The 
intercept between the flat curve and the ideal gas behavior 
provides an accurate estimate of the critical micelle concen- 
tration. 



trends are seen at other T or p. 

An indication of the relative stability of the clusters as 
a function of their size is provided by the cluster potential 
energy. Fig.[7phows the average potential energy (E(s)), 
averaged over all clusters of the same size, vs. size s for 
different T at fixed p. One notices that for each size there 
are several distinct arrangements with different energies 
and entropies. On cooling, lower energy clusters become 
preferentially selected and (E(s)) decreases. One also 
notices a plateau between size 10 and 30 followed by an 
additional plateau, at low T, starting from s ~ 40. The 
initial location of the plateaus coincides with the size of 
the mostly represented clusters. This can be understood 
by considering that the free energy of the system has 
both energetic and entropic contributions. When the po- 
tential energy per particle does not change with the size 
of the cluster, then it becomes convenient to the system 
to favor the formation of the smallest possible cluster 
with the same energy, since it will be possible in this way 
to maximize the total number of clusters and hence the 
translational component of the entropy. 

Fig. [7] also shows the lowest energy configuration 
Emin(s) found for each cluster size, independently from 
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FIG. 6: Cluster size distributions at p = 0.001 for several 
T values (a) and at T = 0.25 for several densities (b). The 
distributions have been normalized so that J2 S s N(s) = 1. 



temperature and density of the simulation. This curve 
provides an estimate for the cluster ground state and 
clearly show that the lowest-energy structures are the 
vesicles, i.e. clusters composed between 40 and 60 parti- 
cles. 



We have examined the connectivity properties of the 
different state points by evaluating the presence of clus- 
ters spanning the entire simulation box. When more than 
50 per cent of the configurations are characterized by 
spanning clusters, the state point is classified as perco- 
lating. Fig. [4] shows the location of the percolating state 
points in the phase diagram. The cluster gas phase is 
never percolating, while the liquid states always are, so 
that percolation properties can also be used to distin- 
guish between the two phases. We also confirm that 
the critical point is located inside the percolation re- 
gion, confirming once more that a pre-requisite for criti- 
cal phenomena is the existence of a percolating network 
of interactions [48 . 
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C. Structure Factor 

It is particularly relevant to look at the structure fac- 
tor in this system both in the gas phase at low T, when 
micelles and vesicles become prominent, as well as in 
the liquid phase. Fig. [8} (a) shows the evolution of the 
structure factor on increasing density at T = 0.25, across 
the gas-liquid transition. At this low T, as discussed in 
Sec. |IVB[ vesicles are the dominant clusters. The struc- 
ture factor indeed evolves from the one characteristic of 
an ideal gas of spherical vesicles (and indeed the inset 
shows that S(q) is well represented by the form factor 
P(q) of a sphere of radius R) 



p(q) 



3(sin(#i?) — qRcos(qR)) 



(4) 



to the one of interacting spheres, in the micelle-rich gas 
phase. Indeed, on increasing p, oscillations at qa ~ 1.2 
arise, which correspond to distances comparable to the 
vesicle size. Beyond density p = 0.6 the gas condenses 
into the liquid phase, where only a very weak pre-peak 
around qa ~ 3 is found. 

Fig. [8]-(b) shows the evolution of the structure factor 
on cooling along the p = 0.4 isochore. Compared to 
simple liquids, one observes a non negligible scattering 
at small q, which progressively increases on approach- 
ing the phase separation. These are the standard critical 
fluctuations which are expected to diverge on approach- 
ing a spinodal line. While in simple liquids, below the 
spinodal temperature the system phase separate in a gas 
coexisting with a liquid phase, here, the peculiar shape 
of the gas-liquid coexistence line (Fig. [l]) opens up new 
stable states, composed by interacting vesicles and S(q) 
becomes peaked at the vesicle-vesicle distance. 
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FIG. 8: Structure factor at (a) T=0.25 for several densities 
and at (b) p — 0.4 for several T. The inset in (a) shows the fit 
with the form factor of a sphere (Eq.El. The best-fit value for 
the radius is R — 2.26a. The above isothermal and isochore 
paths are indicated with dotted lines in Fig|4] 



each other (as in a double-layer), while the peak at fix n 2 
close to 0.9 arises from particles which are forming the 
shells of the vesicle. 

We note that, to quantify the orient ational ordering in 
the system, it is necessary to identify an order parameter 
which does not depend on single-particle properties, since 
both gas and liquid phases are isotropic. We have also in- 
vestigated the rotationally invariant local order [49 indi- 
cators, which have been often exploited to quantify order 
in crystalline solids, liquids and colloidal gels [50, 51 . 
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where qi m (i) is defined as, 



qim(i) 



N, 



1 **< 



(6) 



i=i 



Here A/5, is the set of bonded neighbors of a particle 
i. The unit vector r^- specifies the orientation of the 
bond between particles i and j. In a given coordinate 
frame, the orientation of the unit vector f ^ uniquely de- 
termines the polar and azimuthal angles Oij and <pij . The 
Yim{6ij,4>ij) — Yim(*ij) are the corresponding spherical 
harmonics. We have calculated the distribution of q\ for 
the present model. As shown in Fig. [9|(b), the result- 
ing distributions show systematic differences between gas 
and liquid phase, but it is hard to provide a clear con- 
nection with the degree of order in the system. 



E. Pathways for vesicles formation 



D. Angular Correlations 

To provide evidence of the different ordering of the 
particles in the gas and in the liquid phase, we have cal- 
culated the distribution of relative orientations between 
all pairs of bonded particles. More precisely, we have 
evaluated the distribution of the scalar product fii • A2 
where fix and 112 are the two unit vectors indicating the 
location of the patch center in each particle frame, for all 
bonded pairs. The distribution P(rii • n 2 ) is expected to 
show well denned peaks for an ordered state and to be fiat 
in a completely disordered state. Fig. [9] shows such dis- 
tribution for several p values at low T. The orientational 
differences between the ordered gas phase and the disor- 
dered liquid phase appear very clearly. On increasing p, 
P(rii • 112) evolves from a highly structured function with 
peaks close to |ni • rx2 1 > 0.5 in the vesicle-rich gas phase, 
to a more uniform distribution for the liquid phase state 
points, going back to an evenly more structured function 
when the lamellar phase is entered. The peak at fix • fi2 
close to —1 is created by the pairs of particle which face 



Despite Monte Carlo simulations do not allow for a pre- 
cise definition of time, it is interesting to analyze how par- 
ticles rearrange themselves into large aggregates. More 
specifically, the analysis of the sequence of MC configura- 
tions may help understanding how vesicles are generated. 
We have examined graphically several processes of forma- 
tion of a vesicle starting from a monodisperse solution 
of particles and one of these processes is documented in 
Fig. 10 In all cases, the system starts by forming small 



clusters that grow by incorporating isolated monomers 
or colliding to similar clusters, to form micelles (Fig. [To} 
(a-c)). The transition to vesicle requires the interaction 
between three micelles (d-e), with the formation of a long 
cylindrical micelle (f) which then self-restructures itself 
into the more energetic vesicles configuration (g). 



V. CONCLUSIONS 

This article reports a numerical study of a simple po- 
tential for Janus-like particles, i.e. colloidal spherical 
particles whose surface is divided evenly into two areas 



0.05 





FIG. 9: (a) Distribution of the scalar product between the 
unit vectors of all bonded pairs of particles. The gas-like 
state points are characterized by peaked distributions, reveal- 
ing the presence of orient at ional order. The vesicle configura- 
tions significantly contribute to the region ni • 112 ~ 1, while 
the micelle configurations contribute to ni 112 ~ — 1. (b) 
Distribution of the rotational invariant qi in the gas and in 
the liquid phase. 



of different chemical composition. This study is moti- 
vated by the ongoing effort in the direction of synthesiz- 
ing optimal Janus colloids [9} [TTHT31 [52] . The possibility 
of creating particles whose surface has different behavior 
on the two hemispheres significantly enlarges the rich- 
ness of the resulting collective behaviors. The phase di- 
agram, the ordered and disordered stable structures, the 
self-assembly properties of optimal clusters are indeed 
expected to finely depend on the chemistry and physical 
properties of the particle surface. For this reason, the 
assembly behavior of Janus particles is receiving a con- 
siderable attention even from a theoretical and numerical 
point of view pIH [20j [53H55] • With proper choices of the 
chemico-physical surface properties, Janus particles can 
provide the most elementary and geometrically simple 
example of a surfactant particle [56], in which solvophilic 
and solvophobic areas reside on different part of the sur- 
face of the same particle. 

Modeling of the phase behavior of these particles can 



be performed at different levels of realism. In this article 
we have chosen to implement the highest coarse-graining 
procedure, by considering the two sides of the particle as 
repulsive and attractive, and modeling each of them with 
the simplest corresponding potential, a hard-core plus a 
square-well. Despite this strong simplification, the phase 
diagram of the system displays a very rich behavior, with 
a colloidal-poor (gas) colloidal-rich (liquid) de-mixing re- 
gion, which is progressively suppressed by the insurgence 
of micelles, thus providing a model where micelle forma- 
tion and phase- separation are simultaneously observed. 
The study of this model shows that the suppression of the 
phase separation is driven by the possibility of building 
low energy clusters which are shielded by the presence of 
an external hard-sphere surface, diminishing the driving 
force for forming large-size aggregates. As a result, at low 
T, the system organizes into a dispersion of orientational 
ordered micelles and vesicles, essentially interacting via 
excluded volume only. 

Another advantage of studying the one-patch Kern- 
Frenkel potential is that the parameters of the potential 
can be tuned continuously from the square-well to the 
Janus potential, simply by changing the bonding angle, 
that is the coverage %. This has made possible to follow 
the evolution of the gas-liquid coexistence curve with x> 
thus including the Janus case which displays the inher- 
intly interesting simultaneous presence of a critical point 
and of a micelle formation. It would be interesting in 
the future to explore more in details the transition from 
the standard gas-liquid coexistence to the more exotic 
Janus case exploring a more refined grid of x values. It 
is of course possible that the stabilization of micelle and 
vesicles can take place already at x — 0-6 but in a low T 
region where numerical simulations can not be performed 
at the present time. 

Several questions are still open including the role of the 
interaction range. Experimentally realized Janus parti- 
cles are characterized by interaction ranges of the order 
of 0.1a or smaller. The range width is an important vari- 
able, since it controls both the internal flexibility of the 
aggregates (and hence the cluster entropy) and the ge- 
ometries of the lowest-energy aggregates (cluster energy) . 
In this respect, it can not be a priori foreseen if micelle 
and vesicles will remain the most stable clusters on fur- 
ther decreasing the range. We are currently repeating 
our study for smaller ranges, but this requires a major 
computational effort, since clustering shifts to lower tem- 
perature, making studies of the phase diagram (in par- 
ticular Gibbs ensemble calculations) extremely hard with 
present time numerical resources. The sensitivity of the 
cluster shape to the cluster range would be an impor- 
tant observation and could teach us a lot about how to 
control the shape of the assembly experimentally, and, 
perhaps, help us understanding geometric arrangements 
which are found in protein aggregates. In this respect it 
is worth noticing that micelles have been detected in the 
most recent experimental work of Granick[ll together 
with Bernal-spiral clusters, which have not been observed 
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FIG. 10: Snapshot from a simulation at p — 0.005 and T=0.25 for several MC steps, indicated in the labels. The initial 
configuration, composed by isolated monomers, quickly evolves to form small micelles (a-c). The final vesicle (g) arises from 
a collision between three distinct micelles (d-e) which form an elongated transient tubular cluster (f). Each picture has a side 
length of 9a. 



with the present model, perhaps due to differences in the 
interaction ranges. 
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